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EA ble investment decisions and neglects flexible designs. A more comprehensive approach should address 
Mixed-integer nonlinear programming flexibility in designs. i i . . 
Hydropower This paper proposes an advanced RO methodology for RE generation projects planning, and illustrates 
the methodology using variations of a hydropower case study. The fundamental differences between 
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Nomenclature 


a head (m) 

b scenario probability 
B benefits function 
c costs factor 

C costs function 
Ct costs (109 $) 

d discount rate (%) 


EC,ec economic constraints 

f power plant efficiency 

fe fixed costs (10° $) 

h hours in a year (h) 

H generation capacity (MW) 

K amount of scenarios used for the simulations 
(0) opportunity loss 

p price of electricity ($/kWh) 

P power generated (MWh) 

pf power factor 

pwc present worth factor for costs 

pwr present worth factor for revenue 

R investment timing decisions matrix € [01] 
Rv revenue (10? $) 

TC, tc technical constraints 

To total opportunity loss 


U value according to an evaluation criterion 
V storage capacity (10° m3) 
vcv reservoir’s variable cost (10? $/m?) 
vch power plant’s variable cost (109 $/MWh) 
w additional water inflow (m?/s) 

WwW evaluation criterion weight value 

X water outflow (m? /s) 

y seconds in a year (106 s) 

Y design parameter 


Subscripts 

i index of time 

j index of site 

k index of scenario 

l index of time (Auxiliary) 

m index of evaluation criteria 
n index of investment scheme 
r index of real option 


Superscripts 
Max maximum possible value 
Best best value (higher or lower) according to an evalu- 


ation criterion 


[2], emission targets under the Kyoto protocol [3], renewable sup- 
port schemes such as feed-in tariffs and investment grants [4], as 
well as upgrades in RE technologies achieved through research and 
development [5,6]. However, REs tend to be less competitive than 
other generation technologies, especially in the absence of ade- 
quate support policies [2,7]. This suggests that additional economic 
drivers might be necessary to promote investment in RE generation 
projects. 

A potential and unexploited economic driver for RE generation 
projects is the value of flexibility. In this paper, flexibility refers to 
the capability of managers to modify projects according to the evo- 
lution of uncertainty. This flexibility to adjust projects can enhance 
the projects worth [8], yet itis typically disregarded in the planning 
process of RE generation projects. 


Generation projects, including RE projects, are typically assessed 
with Discounted Cash Flows (DCF) methods [9]. DCF methods 
neglect flexibility, which can lead to project undervaluation. On 
the other hand, RO methods take account of flexibility, and liter- 
ature suggests that RO methods are more suitable for RE project 
assessment [10,11]. Flexibility value captured by RO methods can 
enhance the value of RE generation projects. However, RO liter- 
ature regarding RE generation projects is limited and the theory 
requires further development to take advantage of more flexibility 
value within RE generation projects. 

Flexibility exists whenever projects are influenced by uncer- 
tainty, and means to adjust them (flexibility sources) are available. 
Typical uncertainty sources in RE generation projects are the price 
of electricity, RE source, and technology [12]; whereas flexibility 
sources addressed by RO methods are (i) flexible investment timing 
and (ii) flexible design [13]. Based on these flexibility sources, the 
following classification is proposed for RO planning methodologies: 


e Typical RO: RO methodologies addressing flexibility sources inde- 
pendently. 

e Advanced RO: RO methodologies addressing both flexibility 
sources simultaneously or alternately. 


In other words, typical RO methodologies identify and assess 
flexible investment timing regardless of the projects’ characteris- 
tics (treating projects as black boxes), or flexible projects’ designs 
regardless of possible investment timing. On the other hand, 
advanced RO methodologies identify flexible projects’ designs 
based on possible investment timings and vice versa. 

In this paper, an advanced RO methodology for RE generation 
planning is proposed. The methodology is illustrated and compared 
with other available planning tools using variations of a hypothet- 
ical hydropower case study. 


2. RO theory 


RO theory is the extension of financial options theory for real 
asset assessment. A real option can be defined as the right, without 
obligations, to make an investment decision concerning real assets 
(i.e. defer, build, abandon, alter, switch, etc.). This flexibility can 
enhance the value of projects [14]. 

The evolution of financial options theory, and thus creation of RO 
theory, began in the early 70s due to growing dissatisfaction with 
DCF methods [15], and significant breakthroughs in options pricing 
brought about by the development of the Black-Scholes equation 
[16]. In the late 70s, the term “real option” was introduced, and the 
RO theory was recognized as a new research area [17]. Introduction 
of RO theory triggered research addressing assessment of projects 
in different areas, including RE generation projects. 

Nowadays, relevant RO literature regarding RE and other 
generation technologies is still scarce, and mainly consists of stud- 
ies focused on flexible investment timing that neglect projects’ 
designs. Correia et al. [18] analyze flexible investment timing for 
generation projects. The characteristics of the project are only 
reflected in costs. Chung-Hsiao and Min [19] study the effect of mar- 
ket interdependency on the flexible investment timing of generic 
generation projects. Botterud et al. [20] present a study involving 
different generation technologies. However, the technologies are 
never specified. 

The same trend can be observed in RO literature regarding other 
types of projects, as can be seen in [21-23]. This shows that, cur- 
rent RO literature commonly disregards flexible designs. A reason 
for this tendency can be the underlying difficulties associated with 
assessment of flexible designs. Such assessments require under- 
standing of projects’ technologies, and resulting RO tend to be 
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path-dependent and highly interdependent [24]. Regardless, a few 
examples of literature addressing flexible designs can be found. 
Examples of such literature include [25,26]; where the RO analyses 
are greatly simplified by defining flexible designs empirically. 

The next step in RO theory evolution would involve methods 
for flexible design formulation, and the combination of flexible 
designs and investment timing. A first attempt to take such a step 
was the introduction of a RO approach that, for simplification, will 
be referred to as Tao Wang’s Methodology (TWM) [27]. TWM can 
formulate flexible designs, but only addresses flexible investment 
timings and designs independently. A more advanced approach 
would need to assess both flexibility sources simultaneously or 
iteratively. 

Advanced RO approaches have potential to explore flexibil- 
ity value more than current RO methodologies. Accordingly, they 
could enhance significantly the value of RE generation projects. An 
advanced RO planning methodology is presented in the following 
section. 


3. Methodology 


The proposed RO planning methodology is divided into two 
stages. First, a mathematical programming model, based on a path- 
dependant scenario tree, is used to identify optimal investment 
timings and designs. The model’s output is an investment scheme 
based on the decision nodes considered for the scenario tree. Later, 
Monte Carlo simulation [28] is used to compute the expected 
value of the project when implementing the investment scheme. If 
several investment schemes are proposed, minimax regret (oppor- 
tunity loss) [29] is used to differentiate them. 

The optimization model’s objective function involves profits 
maximization throughout the path-dependant scenario tree. The 
model’s variables are investment timing, and design parameters, 
as can be seen in (1): 


Maximize: S be) 9 [Bir(Yirs Reir) Cir Ykr Rear) (1) 
k i r 


The model’s constraints are based on economic and technical 
limitations of the project, as shown in (2) and (3): 


TC(Y,, Ry) = tc (2) 
EC(Yx, Ry) = ec (3) 


Simulation requires investment schemes and stochastic mod- 
els of relevant uncertainty sources such as, for example, electricity 
prices and environmental policies. The outputs are expected imple- 
mentation values, and a comparative study (based on opportunity 
loss) of the different investment schemes analyzed. Opportunity 
loss calculation is based on selected decision criteria, for example 
net present value and internal rate of return, and a weighted sum 
of objective functions method [30]. Investment scheme selection 
is based on minimum opportunity loss. These computations are 
shown in (4) and (5): 


Onn = 2 >. [URE — Usman (4) 
k 


Ton = S Win Onn (5) 
m 


As can be seen, the models require specific economic and tech- 
nical constraints related to the projects. The formulation of such 
constraints will be illustrated in the case study presented in the 
next section. 


Hydropower plant 


Site 1 
pEr dropower plant 
VA with storage 
— River 

Site 2 == Tunnel 


— Main water flow 
--» Additional water flow 


Fig. 1. River basin. 


4. Case study 
4.1. Description of the case study 


A Generation company can build hydropower plants, and its cor- 
responding dams or tunnels, in three locations along a river basin 
as shown in Fig. 1. Each location has different characteristics which 
affect the design limits and costs of possible hydro power plants. 
The company will assess investments based on expected profits. 

The project’s budget allows construction of at most one power 
plant every 10 years. Construction cannot be delayed more than 
20 years. Operation lifetime is assumed to be 60 years. Relevant 
sources of uncertainty are the price of electricity (modelled with 
Geometric Brownian Motion), and water flows (modelled with log- 
normal distributions). For further details on the models refer to 
[27]. 


4.2. Simplification of the case study 


A simplification of the case study is used to illustrate the main 
differences between traditional DCF, typical and advanced RO 
methodologies. The simplifications are the following: 


e The river basin is simplified as is shown in Fig. 2. 

e All sites’ design limits and storage capacity are assumed to be the 
same. 

e Water flow stochasticity is neglected. 

e An operation lifetime of 30 years is assumed. 

e Operation begins at full capacity five years after construction. 

e Costs for site two and three are assumed 0.5 and 1.0% more expen- 
sive than costs for site one. 

e All parameters used for this case study are given in Table 1. 


The optimization model adjusts the design variable (H) and 
investment timings (R) to maximize profits. This is shown in 


(6)-(8): 


Maximize : X “by (0.001 Ruy — 1000Ct;) (6) 
k 


[Hydropower plant 
= River 
— Water flow 


Fig. 2. Simplified river basin. 
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Table 1 
Parameters for the simplified case study. 
Parameter Value Units 
a 200 m 
b [0.139 0.233 0.233 0.395] 
d 8.6 % 
p [36.24 36.24 36.24 36.24 45.18 45.18 29.11 29.11 56.29 36.24 36.24 23.31] $/kWh x 10-3 
f 0.7 
fe 1.0161 $ x 109 
h 8760 h 
HMax 3200 MW 
pf 0.35 
pwc [1 0.438 0.192] 
pwr [2.8256 0.00 0.00 3.1088 1.2383 0.00 1.4465 1.9963 1.4175] 
y 31.536 sx 10° 
vcv 0.00008070 $/m? x 10? 
vch 0.00010632 $/MWh x 109 
yMax 12,500 m? x 10% 
w [59000] m3/s 


i 
Ruy = SOS PijaPik > PWT AR Lk (7) 
ij I=1 


Ct, = X pwa Ringe + ucvV + vchHhę j) (8) 
i j 


At most, one hydropower plant can be built every time period 
(10 years), as shown in (9): 


X Rij <1 Vik (9) 
i 


Generation (H) and storage (V) capacities have maximum limits 
as shown in (10) and (11): 


V < vMax (10) 
Hj < HV% Vik (11) 


Power generation is a function of generation capacity, water 
flows, dam shape, and storage as shown in (12)-(14): 


Pi jk <2.73fVXij na? Vij,k (12) 
Pi jx < pfhHjk VYij,k (13) 
V > 0.140? (14) 


Construction and operation of hydropower plants affect water 
flow as shown in (15) and (16): 


Xiok =9 Vik (15) 


i 
0.25VS Rij < (Wj +Xij-1,k ~ Xij) Y VYij,k (16) 
l=1 
Timing and design options are formulated based on a binomial 


path-dependent scenario tree addressing three time periods, as 
shown in (17)-(19): 


Ry Aj = R1 j,2Hj,2 = R1 j,3Hj,3 = R1 j4Hj4 Vj (17) 
R2 j, 1Hj1 = Roj2Hj2 Yj (18) 
R2j,3Hj,3 = R2 j,4Hj4 Yj (19) 


The model is a mixed integer non linear programming model. 
Such models can be solved with non linear optimization software, 
such as GAMS®, or a combination of linear optimization software 
and linearisation tools, such as Xpress-IVE® and separable pro- 
gramming techniques [31]. Additionally, due to the size of the case 
study, the model can be treated as a linear optimization prob- 
lem. This can be achieved by conducting an exhaustive search [32] 


for all possible investment decisions (R) (2715 combinations), and 
solving them individually with linear optimization software. The 
exhaustive search can be conducted using MATLAB®, whereas the 
optimization can be performed with Xpress-IVE®. 


5. Evaluation of the simplified case study 
5.1. Traditional DCF analysis simulation 


DCF analyses assume irreversible investments without delay 
possibilities (now or never decisions) [14]. To simulate a DCF 
approach, a greedy optimization [32] based on the net present value 
(NPV) calculation is used. First, optimal design parameters for all 
hydropower plants are determined for the first time period using 
the optimization model. Designs are shown in Table 2. Afterwards, 
investment timing is determined with the greedy optimization. The 
resulting investment scheme can be seen in Fig. 3. 

Based on this methodology, the simplified hydropower project 
expected NPV is 70$ x 10°. 


Table 2 
Fixed design parameters for the simplified case study. 
Site H (MW) 
1 2200 
2200 
3 2200 
56 $/MWh 


aa As 
45 $/MWh Site 3: 2200 MW 


Site 2: 2200 MW 
36 $/MWh 
Site 1: 2200 MW 


36 $S/MWh 


29 $/MWh 


23 $/MWh 


0 10 20 
YEAR 


Fig. 3. DCF method's investment scheme for the simple case study. 
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56 $/MWh 
Sa. a 7 
45 $/MWh Site 2: 2200 MW 


Site 1: 2200 MW 


36 $/MWh 


Site 1: 2200 MW 


36 $/MWh 


29 $/MWh 


23 $/MWh 


0 10 20 
YEAR 


Fig. 4. Typical RO method’s investment scheme for the simple case study. 


5.2. Typical RO simulation 


To simulate a typical RO analysis, design flexibility is neglected. 
A typical RO analysis considers flexible investment decisions by 
assessing the option to delay investments and gain information 
about the evolution of uncertainty. To simulate this approach, 
design parameters are fixed to those given earlier in Table 2, 
and optimal investment timing is determined with the opti- 
mization model. The resulting investment scheme is shown in 
Fig. 4. 

The results show that it is convenient to delay investment for at 
least one time period. This flexible investment decision enhances 
expected NPV from 70 to 97$ x 10°. 


5.3. Advanced RO simulation 


Advanced RO analyses exploit the value of flexible investment 
timing and designs. This considers the value of delaying investment 
and adjusting design parameters based on selected investment 
timing. This approach is simulated using the optimization model 
presented in the previous section. The resulting investment scheme 
is shown in Fig. 5. 

The advanced RO methodology also proposed delaying invest- 
ment until the second period. In addition, it suggested a lower 
generation capacity for the possible hydropower plant located 
on site two (from 2200 MW to 2000 MW). This flexibility, in the 


56 $/MWh 
ER 7 
45 $/MWh Site 2: 2000 MW 


Site 1: 2200 MW 


36 $/MWh 
Site 1: 2200 MW 


36 $/MWh 


29 $/MWh 


23 $/MWh 


0 10 20 
YEAR 


Fig. 5. Advanced RO method's investment scheme for the simple case study. 


56 S/MWh 


45 $/MWh 
3600 MW 


Site 2396 310° 


36 $/MWh Wait 


Site 3: 1700 MW 
Tunnel 


36 S/MWh 
Wait 


29 $/MWh Wait 


23 SIMWh 


0 10 20 
YEAR 


Fig. 6. DCF method’s investment scheme for the case study. 


project’s design and timing, increases the expected NPV from 97 to 
98$ x 10°, 

It can be argued that this study does not show significant 
economic advantages of advanced RO methodologies over typical 
RO approaches. Nevertheless, it is a simple and replicable study 
that illustrates the fundamental differences between the three 
approaches. The economic advantages of advanced RO methodolo- 
gies will be addressed in the following section. 

To summarize this section, traditional DCF methodologies 
focus on selecting the best investments for the time being; 
typical RO approaches enhance the value of projects by either 
addressing flexible investment timing, as shown in the case 
study, or flexible designs; and advanced RO methodologies fur- 
ther enhance the value of projects by considering flexible designs 
and timing simultaneously, as shown in the case study, or alter- 
nately. 


6. Evaluation of the complete case study 
6.1. Investment scheme formulation 


In this section, the complete case study is assessed with 
available approaches representing DCF, typical and advanced RO 
planning methodologies. DCF methodologies are represented by 
a greedy optimization based on the calculation of NPV. Typi- 
cal RO methodologies are represented Tao Wang’s Methodology 
(TWM) [27], which is the most comprehensive typical RO plan- 
ning methodology found in the literature survey. Advanced RO 
approaches are represented by the methodology proposed in this 
paper. 

In order to define project designs for the DCF analysis, part of 
TWM is used. The investment scheme formulated with the DCF 
methodology is shown in Fig. 6; whereas the designs used for the 
analysis are presented in Table 3. 

First, TWM formulates several design alternatives using an opti- 
mization model based on possible values of underlying sources of 
uncertainty. The designs are given in Table 4. A second optimiza- 


Table 3 
Fixed design parameters for the case study. 


Site Generation (MW) Storage (m? x 10°) 
1 3200 12,500 
2 3600 9600 
3 1700 0 


2092 E.A. Martinez-Cesefia, J. Mutale / Renewable and Sustainable Energy Reviews 15 (2011) 2087-2094 


Table 4 
Design parameters formulated with TWM. 


Design parameter Site First alternative Second alternative Third alternative 
Generation (MW) 1 1723 1946 1966 
2 3600 3600 3600 
3 1700 1700 1700 
Storage (m? x 10°) 1 9593 12,242 12,500 
9600 9600 9600 
3 0 (tunnel) 0 (tunnel) 0 (tunnel) 

56 $/MWh 56 $/MWh 


1966 MW 
12.5 m3103 


36 $/MWh 


23 S/MWh 


0 10 20 
YEAR 


45 $/MWh Site 1: 


3600 MW 
9.6 m310° 


Site 2: 


36 S/MWh 
Site 3: 1700 MW 


Tunnel 


29 $/MWh 


Fig. 7. TWM'’s investment scheme for the case study. 


tion model is then used to select a design alternative and calculate 
its optimal investment timing. The resulting investment scheme is 
shown in Fig. 7. 

Table 5 shows different design alternatives for site one; whereas 
the design alternative remains unchanged for other sites. This indi- 
cates that flexible designs were only identified for site one. 

The proposed advanced RO methodology adjusts design param- 
eters and searches for optimal investment timing in a single 
process. The investment scheme formulated with this methodology 
is shown in Fig. 8. 


6.2. Investment scheme evaluation 


Investment schemes are evaluated using Monte Carlo simula- 
tion (5000 scenarios) and opportunity loss criteria. The simulation 
requires investment schemes (formulated in the previous sub- 
section), and stochastic models for relevant sources of uncertainty 
namely electricity price and water flows. 

The electricity price is modelled annually using Geometric 
Brownian Motion and considering an initial value of 36.24 $/MWh, 
a drift rate of 0.33% and a volatility of 6.96%. The water flows are 
modelled seasonally (every six months) with lognormal distribu- 
tions. The main river's mean water flow is 374 and 283 m?3/s, witha 
standard deviation of 87.7 and 45.4% respectively. The mean water 
flow of the additional river is 389 and 154 m?/s, with a standard 
deviation of 89.8 and 9.7% respectively. 


45 $/MWh 


3600 MW 
te 2: 
Site 2:3 45 m3109 


36 $/MWh 


Site 3: 1700 MW 
Tunnel 


3600 MW 


Site 259 6 m10? 


29 $/MWh 


23 $5/MWh 


0 10 20 
YEAR 


Fig. 8. Advanced RO’ investment scheme for the case study. 


Opportunity loss is based on NPV. However, for illustrative pur- 
poses, the net value (NV) will also be employed. The simulation 
results are shown in Tables 5 and 6. 

As can be seen in Table 5, the advanced RO investment scheme 
offers the highest expected value, followed by investments deter- 
mined using typical RO and DCF methodologies. 

It can be seen, in Table 6, that the opportunity loss of the 
advanced RO investment scheme is greater than zero, but sig- 
nificantly the lowest. This suggests that although the investment 
scheme is not the optimal in all possible scenarios, it is the opti- 
mal or close to the optimal investment scheme in most scenarios. 
It can also be seen that NV rather than NPV loss of opportunity is 
more favourable for RO investment schemes. To explore this fact, 
the NV and NPV probability distribution functions are shown in 
Figs. 9 and 10. 

The curves behaviour is more similar for NPV than for NV, this is 
caused by the time value of money. In this case study, RO for flex- 
ible investment timing balance timing benefits (additional profits) 
against timing costs (time value of money). 

Therefore, presented results, besides being dependent on the 
stochastic models selected, can be significantly dependent on the 
discount rate. 


6.3. Sensitivity analysis 


To address the robustness of the presented results, a sensitivity 
analysis is conducted for the Monte Carlo simulation. The analysis 


Table 5 Table 6 

Mean expected NPV and NV for all investment schemes. Mean expected opportunity loss for all investment schemes. 
Investment scheme NPV ($ x 10°) NV ($ x 10°) Investment scheme ($ x 10°) NPV regret ($ x 10°) NV regret 
DCF 1159 19,595 DCF 34 8805 
Typical RO 1160 22,556 Typical RO 325 5844 
Advanced RO 1180 28,252 Advanced RO 13.3 148 
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—+— DCF 


Typical RO 
—— Advanced RO 


Probability (%) 


NV ($ 109) 


Fig. 9. Net value probability distribution function. 


—+— DCF 


Typical RO 
—— Advanced RO 


Probability (%) 


5 10 15 20 25 
NPV ($ 108) 


Fig. 10. Net present value probability distribution function. 


Table 7 
Probability of being the best or worst investment scheme. 


Investment scheme NPV best(%) NV best (%) NPV worst (%) NV worst (%) 


DCF 0.0 0.0 59.2 100 
Typical RO 0.0 0.0 40.7 0.0 
Advanced RO 100 100 0.0 0.0 


considers 10% variations in discount rate, electricity price model’s 
drift and volatility, and average water flows. The results are shown 
in Table 7. 

It can be seen that, within the scope of the sensitivity analysis, 
the advanced RO investment scheme is the best alternative. 


7. Conclusion 


This paper has discussed the importance of RE generation 
projects, as well as the necessity for additional economic drivers 
for investments in RE. RO theory is explored as a possible eco- 
nomic driver for RE projects. It is concluded that RO theory requires 
further development to significantly enhance RE projects. Conse- 
quently, advanced RO methodologies are suggested to provide this 
improvement. 

An advanced RO methodology is introduced and illustrated with 
variations of a hydropower case study. 


A simplification of the case study is used to illustrate the 
theoretical difference between DCF, typical and advanced RO 
methodologies. In the absence of either flexible designs or flexible 
investment timing, an advanced RO methodology behaves as a typ- 
ical RO approach. Moreover, in the absence of all flexibility, a typical 
RO approach behaves as a DCF technique. This suggests that the 
lowest value an advanced RO methodology would provide would 
be the optimal value according to a traditional DCF methodology. 

The complete case study allowed a comparison of the proposed 
advanced RO methodology against other planning methodologies. 
Results showed higher expected profits for projects planed with 
the advanced RO methodology. 
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